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Abstract 

We analyze the effects of thermal conduction in a relativistic fluid just 
after its departure from spherical symmetry, on a time scale of the order of 
relaxation time. Using first order perturbation theory, it is shown that, as 
in spherical systems, at a critical point the effective inertial mass density 
of a fluid element vanishes and becomes negative beyond that point. The 
impact of this effect on the reliability of causality conditions is discussed. 



1 Introduction 

The behaviour of dissipative systems at the very moment when they depart 
from hydrostatic equilibrium has been recently studied ||, || . 

As result of these works, it appears that a parameter formed by a specific 
combination of relaxation time, temperature, proper energy density and pres- 
sure, may critically affect the evolution of the object. 

*On leave from Departamento de Fisica, Facultad de Ciencias, Universidad Central de 
Venezuela, Caracas, Venezuela and Centro de Astrofi'sica Teorica, Merida, Venezuela. 
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Specifically, by means of first order perturbation theory, it was shown that 
in the equation of motion of any fluid element, the inertial mass density term 
is multiplied by a factor which vanish for a given value of the parameter a 
defined below in equation (50) (critical point) and changes of sign beyond that 
value. Nevertheless, even though self-gravitating systems seem to become more 
unstable as the above mentioned parameter grows an exact description of 
the evolution B docs not indicate any anomalous behaviour of the system at 
or beyond the critical point. Thus, first order perturbation theory seems to be 
not reliable close or beyond the critical point. 

Although causality conditions prevent, in some cases (pure shear or bulk 
viscosity reaching the critical point. This may not be the case in the most 
general situation (heat conduction plus viscosity || or even in the non-viscous 
case Q). It is worth mentioning that causality conditions were found using a 
first order perturbation method [[|, thus, in some cases these conditions must 
be revised. 

However in all these works it has been assumed that spherical symmetry 
is preserved. It is therefore pertinent to ask if the appearance of the critical 
point is closely related to that kind of symmetry (spherical) or if it represents 
a general feature of dissipative systems. 

It is the purpose of this work to provide an answer to the question above by 
considering a system which, although initially spherically symmetric, is submit- 
ted to perturbations deviating it from spherical symmetry. 

By examining the equation of motion of an arbitrary fluid element along the 
"meridional" direction, it is shown that a critical point also appears in this case, 
indicating thereby the independence of this effect from the spherical symmetry. 

The paper is organized as follows. In the next section the field equations and 
conventions are presented. In section 3 we briefly present the equation for the 
heat conduction and write down explicitly the ^-component of this equation. In 
section 4 the full system of equations is evaluated at the time when the object 
starts to depart from spherical symmetry. Finally, a discussion of results is given 
in the last section. 

In the whole text, a bar over a quantity denote that this one has been 
measured by a Minkowskian observer. If the symbol over a given quantity is a 
tilde, then this one is evaluated after the perturbation. 



2 The field equations and conventions 

We shall consider axially and reflection symmetric distributions of fluid dissi- 
pating energy through a heat flux vector. In null coordinates the metric is given 
by the standard Bondi expression ^ 

ds 2 = -U 2 r 2 e 2 ^)du 2 +2e 2 Pdudr 

+2Ur 2 e 2 ~>du d9 - r 2 (e^dO 2 + e" 2 T sin 2 9 d(j) 2 ) , { ' 
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where V, (3, U and 7 are functions of u, r and 9. 

We number the coordinates x 0,1,2 ' 3 — u, r, 9, <f>, respectively, u is a time-like 
coordinate such that u — constant defines a null surface. In flat space-time 
this surface coincides with the null light cone open to the future, r is a null 
coordinate (g rr = 0) and 9 and <p are two angle coordinates (see for details). 

In order to guarantee the regularity of the metric within the fluid dis- 
tribuition at r = 0, we impose the conditions Js) 

V = r + 0{r 3 )- /3 = 0(r 2 ); 

U = 0{r); 7 = 0(r 2 ). [> 

In this case, the energy-momentum tensor takes the usual form 



Tfj, u = {p + - pg^f + q^Vv + q v v^ (3) 

where p, p and q M denote proper energy density, pressure and the heat flow 
vector, respectively, and the four-velocity satifies the conditions 



^% = 1 



(4) 



Since reflection symmetry is preserved, we have 



v 3 = 



Q3 = 



(5) 



We shall now write the field equations according to the scheme presented by 
Isaacson et al. ^ (see also J9|). Thus, there are three main equations which 
only contain derivatives with respect to r and 9. In our case the corresponding 
expressions are given by 



1 2 

2nr [{p + p)vi+ 2gi] v 1 = /? 1 - -r (7,1) 



(0) 



167rr 2 [(p + p) viv 2 + q\v 2 + q 2 Vi] 



r 4 e 2(7-/3) [/l 



2r z 



r 2 (r- 2 (3) 12 -sin- 2 #(7f 



2 7,i7,: 



(7) 



2V^ 



-8 7 rr 2 e 2 ' 3 [p - p + r^e^ 7 ([p + p] (v 2 ) 2 + 2q 2 v 2 )] = 
-e 2 ^- fi \U tl ) 2 - (r 2 sinfl)" 1 (r 4 C/ sin 9), 12 + 2e 2 ^ x 



-l + (sin0) 1 (/3,2sin6») !2 -722-37,2 cot 6>+(/3, 2 ) 2 + 27,2(72-/3,2)! (8) 



where a comma denotes partial derivative. Following the Bondi-Isaacson 
scheme we shall assume that variables p, p, Vi, v 2 , q\ , q 2 and 7 are known on 
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a given initial null hypersurface (u = uq = constant), which we shall call the 
"initial cone" . 

At u = uq perturbations are introduced into the system, forcing it to deviate 
from spherical symmetry. We shall evaluate the system immediately after per- 
turbing it, on a null surface (u = u = constant) such that u — uq is much smaller 
than the typical time required for V\, v 2 , qi, q 2 and 7 to change significantly 
(where "significantly" refers to terms linear in the perturbation, or larger). In 
other words, on the surface u = u, which we shall call the "evaluation cone" , 
the components v% and v 2 of the four- velocity, q\ and q 2 of the heat flow vector 
and the metric function 7 have the same magnitude as on the initial cone (up 
to cuadratic and higher terms in the perturbation parameter). 

The "standard equation" which provides an expression for 7.0 reads 



4r(r 7 ) 01 = [2rj A V ~r 2 (2 j, 2 U + U, 2 -U cot 9)] 1 - 2r 2 sin" 1 6» (71C/ sin6») 



Next, in order to give physical meaning to quantities appearing in (|3|), taking 
into account the axial symmetry of the case considered, we shall develop a 
procedure similar to that used by Bondi JhJ in his study of non-static spherically 
symmetric sources. Thus, let us introduce purelly local Minkowski coordinates 



+ ?_ e 2(7-/3) ([/l) 2 +2e 2(/3- 7 ) x 

(A2) 2 + A22 - 13,2 cot 6 + Air ([p + p] {v 2 ) 2 + 2q 2 v 2 ) 



(9) 



(x ,^ 1 ,^ 2 ,^ 3 ), defined by (IT) 



dx° = Adu + Bdr + CdO 



(10) 



dx 1 = Bdr + Cd8 



(11) 



dx 2 = FdO 



(12) 



dx 3 = Gdcj) 

where cofficients A, B, C, F and G are given by 



(13) 




(14) 




(15) 




(16) 
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F = re 7 



(17) 



G 



so that, 



ds 2 



rf\2 



{dary - {dx i y - (dx 2 y - 



(18) 



(19) 



Now, for the Minkowski observer defined above the fluid moves in such a 
way that the four-velocity vector of the fluid is given by 



B* , = A(l > wi,wa,0). 



(20) 



where u>i and 0J2 are the components of the velocity along the x 1 and x 2 
axes, respectively (as measured by our locally Minkowski observer), and 



A 



(21) 



Thus, a locally Minkowski observer comoving with the fluid is one for which 

U>\ = Ll>2 = 0. 

Now, from @ and ©-(H) we can obtain expressions for the four-velocity 
components in the original Bondi coordinates. In terms of u>i and u>2, they are 



A 



1-wa 1 



.4 



or, for the covariant components 

'V 

vq = 
Vi 



v 2 = A 



A ^e 2 ' 3 - U 2 r 2 e 2 ^ 
\e 2 ? (1 - w{) 

(¥. e 2l3 _ U2 r 2 e 2 7 y 

Ur 2 e 2 ~< (1 - Wi) 



W2 

— 



(Z e 2/3 _ [/2 r 2 e 2 7 
«3 - 



1/2 



(22) 

(23) 
(24) 

(25) 
(26) 



It is worth noticing that the Minkowski observer will say that the fluid is at 
rest when u>i — L02 — 0. Also, observe that uj 2 , as it follows from (|lj), measures 
the velocity of a fluid element along the 9 (meridional) coordinate line. 

Finally, defining q as 



(27) 
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we obtain 



q = q\uJi \l -e 2 P - U 2 r 2 e^ (28) 



r 



qi = -qe 



g 2 = -gre 7 x (30) 
C/re- f - UWA f M±^# - Au^ + -J**, (A - 1) 

93 = (31) 

3 Heat Conduction Equation. 

In the study of star interiors it is usually assumed that the energy flux of radi- 
ation (and thermal conduction) is proportional to the gradient of temperature 
(Maxwell- Fourier law or Eckart-Landau in general relativity). 

However, it is well known that the Maxwell- Fourier law for the radiation flux 
leads to a parabolic equation (diffusion equation) which predicts propagation 
of perturbation with infinite speed (see |^2|- |l^| a nd references therein). This 
simple fact is at the origin of the pathologies || found in the approaches of 
Eckart and Landau [[l6| for relativistic dissipative processes. 

To overcome such difficulties, different relativistic theories with non- vanishing 
relaxation times have been proposed in the past ]T^| - |20[ . The important point 
is that all these theories provide a heat transport equation which is not of 
Maxwell- Fourier type but of Cattaneo type |2l| , leading thereby to a hyperbolic 
equation for the propagation of thermal perturbation. 

Accordingly we shall describe the heat transport by means of the relativistic 
Israel-Stewart equation |I3 , which reads 

+ <f = kP ° P {T ' ~ Ta0) ~ TyaqpaP ~ \ kT2 -p qa (32) 

where k, t, T and a' 3 denote thermal conductivity, thermal relaxation time, 
temperature and the components of the four-acceleration, respectively. Also, 
P a/3 is the projector onto the hypersurface orthogonal to the four velocity v a 
and 

Dq a a . 



G 



For the purpose of this work we only need the components of ( p2| ) containing 
u-derivatives of q 2 and v%. Such terms only appear in the 0-component of (p2|). 
Thus we have 

r q v P» 2 + q 2 = nP a2 (T a - Ta a ) - X -kT 2 (^vP) ^ ? (33) 
or simbolically 

Tg ai viq 2 flP 22 + KTP 22 g 01 v lV2fl = H (34) 

where TL is a combination of terms not containing 52,0 or V2,o- 
Alternatively, we may write from ( |34| ) 

kT 

92,o = v 2 . a +J, (35) 

r 

where Z, as Ti., represents a combination of terms not containing 52.0 or t>2.o- 
We shall now get into the central problem of this work. 

4 Thermal conduction and departure from hy- 
drostatic equilibrium 

Let us now consider a fluid distribution which on the initial cone (u — uq) is 
spherically symmetric. We shall not impose any further restrictions on it, so 
that in principle the system may not be in hydrostatic equilibrium along the 
radial direction at u = u . 

Next, let us assume that the system is perturbed, as a result of which it 
departs from spherical symmetry. 

Since the system is initially spherically symmetric it is clear that on the 
initial cone the fluid is in equilibrium along the 0-line of coordinates (i.e. at 
u = mo, ^2 = <72 = and 1^2,0 = 92,0 = 0). 

For the purpose of this paper (i.e. in order to put in evidence the existence 
of a critial point along a non-radial direction) it will suffice to consider those 
equations containing (72.0 and ^2,0- 

Let us now evaluate our system on the evaluation cone [u = u). 

As mentioned in the precedent section, the only component of ( |32| ) containing 
terms with u-derivatives of q 2 or lj 2 (the 9 -component) leads to eq.(^5|). 

This last equation after perturbation keeps the same form as (^5|) , where all 
the quantities are evaluated on the u = u cone, i.e. 

kT „ ~ 

92,o = V2.0+1- (36) 

T 

Next, instead of dealing with the field equations (H)-(||), it will be more 
useful to consider the "conservation" equations 
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= (37) 
Thus, the it-component of ( p7| ) gives an expression of the form 

[(p + p)v + q ] Vifi + [{p +p)vi+ qi] v ,q + qifiVo + q ofi Vi = J (38) 

where J denotes a combination of terms without it-derivatives of matter 
variables. Observe that d38| ) does not contain terms with u -derivatives of q2 or 
0J2- Also note that all terms containing W2,o, appearing from the it-derivatives 
of A in v\ o, vq q, qi o and Qo,o are multiplied by w 2 and therefore vanish on the 



evaluation cone (see eqs.(f22|)-(31 )) 



Then, pq) may be written alternatively as 



/(wi,o;?,o) = (39) 

where /C, as J , is a combination of terms without u -derivatives of matter 
variables. 



The r-component of (37) yields an expression of the form 



[(P + P) vi + qi] vx.o + qi,oVi = C, (40) 

where C is again a combination of terms without it-derivatives of matter 
variables. 

As in the precedent case, when evaluating ( fief ) at it = u , all terms with u>2.o 
will vanish since they appear multiplied by u>2 , and therefore this later equation 
may be written as 

g(ui,o\q,o) = M (4i) 

M. denoting another combination of terms without it-derivatives of matter 
variables. 

Finally, let us consider the ^-component of ( |37| ) , it leads to an expression of 
the form 



[(p + p)v 2 + q2] vi t o + [(p + p)vi+ qi] v 2 fi + q\flV 2 + ?2,oUi = A/" (42) 
or, using @ and @ 

[{p + p) vi + qi ] v 2fi + qifiVi = AT + h (K,M) (43) 

Since J\f does not contain terms with it-derivatives of matter variables, these 
later terms are absent on the right hand side of (^ ) . 

Before perturbation (on the initial cone) ( f43| ) is just an identity = 0, since 
the system is spherically symmetric. After perturbation (^) keeps the same 
form, with all quantities evaluated at u = u, i.e. 
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[(p + p) vi + qi] i) 2 ,o + 92,0^1 = N + h (JC, M) 
or replacing <f 2 ,o by (p 



p + p 



or, after re-arranging terms 




w 2 ,o+2^i =M + h(K,M) 



yp + p- 

Next, taking into account that 




V2.0 



N+h{K,M) 



Vl 



(44) 



(45) 



(46) 



(47) 



as it follows from (gj) and ( f29| ) imposing W2 = 0, and keeping terms linear 
in the perturbation, we may write (Efl) as 



kT 

T 



p + p q v 2 ,o 



Af +h(K.,M) 



■X, 



Vl 



(48) 



It remains to expand V2,o by using (P5|). One obtains terms with u>2,o, terms 
with u>i t o and terms not containing it-derivatives of matter variables. The latter 
may be included in the right hand of ([is]). Terms with Wi,o ma Y also be trans- 
formed into terms not containig it-derivatives of matter variables by virtue of 
@ and @. Thus, we obtain finally 



Re = -W2,o x 



1 - 



kT 



(p + p- q) 



(49) 



t(p + p- q)_ 

The physical meaning of Rg is quite simple. If the fluid is in equilibrium 
along any meridional line then Rg = 0, thus Rg represents the total force on 
any fluid element in the ^-direction (observe that p^ enters in Rg). 

Therefore ( [49"| ) may be regarded as an equation of the "Newtonian" form 

Force = mass x acceleration 

where by "acceleration" we mean the time-like derivative of 012- 
Thus, if 

kT 

a = — = 1 50 

r{p + p-q) 

the effective inertial mass vanishes and any fluid element is out of hydrostatic 
equilibrium (0*2,0 7^ 0) even though Rg is zero. Note that beyond the critical 
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point (a > 1), the system exhibits an anomalous behaviour: If the total force 
on any fluid element is acting in the positive ^-direction, then the system tends 
to accelerate in the negative ^-direction. 

This is the same critical point already found in the spherically symmetric 
case jl], ^, B|. The only difference with the later case being the appearance 
of q in a. This is due to the fact that in the spherically symmetric case the 
system is initially static (q = 0) or slowly evolving (q is small). However in our 
case, we have considered the more general situation when the system initially 
is spherically symmetric but, although q2 vanishes on the initial cone, q\, and 
therefore q, is not necessarily zero or a small quantity. 

5 Conclusions 

We have described the departure of a fluid distribution from hydrostatic equi- 
librium along the meridional direction, which in its turn implies departure from 
spherical symmetry. It has been shown that a critical point, for which the ef- 
fective inertial mass density of a fluid element vanishes, exists. It has been 
therefore established, that the occurence of such critical point is independent of 
spherical symmetry. 

As in spherically symmetric systems jl], g), condition a — 1 establishes 
the upper limit for which first order perturbation theory can be applied. Thus, 
close to (and beyond) this point, causality conditions obtained from the linear 
approximation should be taken with caution. 
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